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ABSTRACT 

In radio astronomy, reference signals from auxiliary antennae, receiving only the radio frequency interference 
(RFI), can be modified to model the RFI environment at the astronomy receivers. The RFI can then be cancelled 
from the astronomy signal paths. However, astronomers typically only require signal statistics. If the RFI 
statistics are changing slowly, the cancellation can be applied to the signal correlations at a much lower rate 
than required for standard adaptive filters. In this paper we describe five canceller setups; pre- and post- 
correlation cancellers that use one or two reference signals in different ways. The theoretical residual RFI and 
added noise levels are examined, and demonstrated using microwave television RFI at the Australia Telescope 
Compact Array. The RFI is attenuated to below the system noise, a reduction of at least 20 dB. While dual- 
reference cancellers add more reference noise than single-reference cancellers, this noise is zero-mean and 
only adds to the system noise, decreasing the sensitivity. The residual RFI that remains in the output of single- 
reference cancellers (but not dual-reference cancellers) sets a non-zero noise floor that does not act like random 
system noise and may limit the achievable sensitivity. Thus dual-reference cancellers often result in superior 
cancellation. Dual-reference pre-correlation cancellers require a double-canceller setup to be useful and to give 
equivalent results to dual-reference post-correlation cancellers. 

Subject headings: instrumentation: interferometers — methods: data analysis — methods: statistical — tech- 
niques: interferometric 



1. INTRODUCTION 

Radio astronomy is poised to move ahead with a suite of 
new instruments such as the Square Kilometre Array (SKA). 
These instruments will improve sensitivity, resolution, band- 
width, and many other instrument and observational parame- 
ters by more than an order of magnitude. At present, many 
radio astronomy observations are corrupted to at least some 
extent by radio frequency interference (RFI). This interfer- 
ence comes from ground-based communication transmitters, 
satellites, and the observatory equipment itself. With the in- 
crease in sensitivity and frequency coverage of radio astron- 
omy instruments, and with telecommunication signals occu- 
pying more of the spectrum, it is essential to develop ways of 
\ removing or suppressing this RFI. 

In radio astronomy, real-time adaptive filters can be used 
to modify an auxiliary voltage time series (the reference 
signal) so that it cancels RF I from an astr o nomical volt- 
, age time series ([Ba rnbaum & BradlevI II998I : iBoweii 1200 It 
' iMitchell & Robertsonn2005i) . For each voltage sample the fil- 
ters are allowed to slightly vary their internal coefficients in 
order to adapt to any changes taking place in the RFI. If one 
is interested in the power spectrum of the astronomy signal, 
and the filter coefficients stay fairly constant over the inter- 
val in which the power spectrum is estimated, cancellation 
in the post-correlation domain can give better results, since a 
second reference antenna can be used to give complete sup- 
pression of the RFI (only zero-mean random receiver noise is 
add ed to the complex correlations and it will av erage away, 
see iBriggs et"ani2000l: IMitchell & Bowed 1200 Ih . However, 
many applications (particularly in the communications field, 
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but also some astronomy applications), require the recovery of 
the actual symbol stream (i.e., a transmitted sequence of sym- 
bols such as bits or words) from the noisy RF environment, 
which is not retained in post-correlation. Following a sugges- 
tion of Briggs et al. ( 2000), we have devised a modified ap- 
proach that can give improved RFI attenuation in the voltage 
domain. In an attempt to minimise any RFI in an astronom- 
ical voltage series the standard approach is to minimise the 
canceller's output power, which it will be shown means that 
some residual RFI always remains. The modified approach 
that we give here forces the RFI in the output power to zero. 
It results in residual power that is always greater than that of 
the standard approach (output power is no longer minimised), 
but which does not contain RFI. That is, superior RFI cancel- 
lation is obtained, but at the expense of somewhat increased 
thermal noise. 

In the following sections the standard adaptive canceller is 
discussed, and the new approach is introduced. This is fol- 
lowed by an overview of how cancellation can be applied af- 
ter correlations are formed. Residual RFI power and added 
receiver noise are investigated and an example from the Aus- 
tralia Telescope Compact Array is given. 

2. ADAPTIVE CANCELLERS 

2.1. The Model 

In an attempt to remain general, we assume a system of 
one or more radio antennae pointing towards a direction on 
the celestial sphere. Delays are inserted into the signal paths 
so that a wavefront from the chosen direction arrives at the 
output of each antenna simultaneously. The c elestial location 
is known as the p hase tracking centre (Tho mson et al.llT986l : 
iTavlor et al.llI999i) . Given that we will attempt to deal with 
the interference in subsequent parts of the system, separate 
reference antennae are incorporated into the network of re- 
ceivers to observe the RF environment, which typically enters 
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the astronomy signal through the side-lobes of the antennae. 
Signals from astronomy antennae will be referred to as main 
signals, and those from reference antennae as reference sig- 
nals. At each antenna a waveform containing an additive mix- 
ture of all the signals present in the environment is received 
and downconverted to an IF voltage series. This is sampled 
and quantised into a number of digital bits (which we assume 
is sufficient to keep the voltage statistics linear so that quan- 
tisation effects such as intermodulation are negligible, and to 
ensure that the receiver noise and astronomy fluctuations are 
measured even in the presence of strong interference). Each 
main IF voltage series contains three components: a noise 
voltage from the receiving system, n(t); a noise voltage from 
the sky, s{t); and interference, i{t). The sky voltage contains 
the information about the astronomical sources, which are the 
signals of interest. If the interference cannot be removed com- 
pletely, it is desirable to reduce it to less than the final RMS 
noise level (with a negligible - or a t least predictable - effect 
on the astronomv. lBarnbaum & Bra dlev 1998). 

Since an interfering signal is usually incident from a direc- 
tion other than the phase tracking centre, its wavefront will not 
be synchronous at the output of the different antennae. The 
geometric delay, r^, of antenna k - due to the physical sep- 
aration of the receivers - represents the difference in arrival 
time of an interfering wavefront at antenna k and an arbitrary 
reference point (after accounting for the delay needed to track 
the selected field on the celestial sphere). 

As a signal passes through a receiving and processing sys- 
tem it encounters various convolutions and deconvolutions, 
so working in the frequency domain can offer a more intu- 
itive basis for discussion. In the frequency domain the system 
can be represented by complex multiplications and divisions. 
A frequency-dependent coupling term, Gy(i), is used to de- 
scribe the combined complex-valued gain of each receiver 
system and antenna to the interference, including any filter- 
ing (time is included to account for the slow variations im- 
posed as the RFI passes through antenna side-lobes). Using 
upper case characters to denote frequency domain quantities, 
and keeping in mind that this spectral representation comes 
from Fourier transforming each consecutive 1000 or so sam- 
ples of the voltage series, the signal in a quasi-monochromatic 
channel at frequency v is 

K,m = N,..m + G,,^{t)I,e^*--^''> + S,, (1) 
for main antennae and 

K,r = N,^r + G,^.(i)/.e^'^-'-(*\ (2) 

for reference antennae. Note that the phase term due to the ge- 
ometric delay of the interfering signal, 4>u,k{t) = 2TTVTk{t), 
has been kept separate from the G-terms (time is also included 
here to account for changes as either the phase tracking cen- 
tre or RFI transmitter direction change). To remain general 
the G-terms are kept as complex quantities, to allow for any 
effects on the phase that are not due to the geometric delay. 
It is assumed that th ere is only on e interfering signal in a fre- 
quency channel (see lBoweiil200 1' for a discussion of multiple 
interferers), and that there is negligible reference antenna gain 
in the direction of S, i.e., reference antennae do not measure 
signal from astronomical sources. These are important as- 
sumptions, but often quite reasonable (the latter assumption 
is strengthened because the weak astronomy signal enters the 
reference antennae through side-lobes). However, statements 
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Fig. 1. — Schematic of a MKla adaptive canceller. Dashed hnes 
indicate signal paths used to calculate the weights. 

made later in relation to the lack of effect of adaptive can- 
cellers on the astronomy signal rely on the validity of the 
second assumption. If any astronomy signals leak into the 
reference series the achievable sensitivity, dynamic range, as- 
tronomy purity, etc., will all be affected. For example, signal 
leakage could lead to the power of a strong self-calibration 
source changing as the canceller weights change. 

2.2. MKl: Single Reference Adaptive Cancellers 

Adaptive cancellers are usually applie d to broadband IF 
voltage samples in the time d omain (IWidrow & Steamsl 
1 1985; Barnbaum & Bradleyl[T998h . ("Broadband" here sim- 
ply refers to the whole passband, rather than the quasi- 
monochromatic frequency channels.) While the adaptive can- 
celler examples given throughout section |4] have been pro- 
cessed in the time domain, the following analysis is carried 
out in the frequency domain ( see Widrow & Stearns 19851: 
iBarnba um & Bradley 19981; and lBowe3l200ll for descriptions 
of time domain implementation). 

The aim of adaptive cancellers in interference mitigation is 
to find the set of filter weights, W, which scale and phase shift 
each reference antenna frequency channel so that they best 
approximate the RFI in a main astronomy spectrum (in the 
time domain delays are inserted into the signal paths so that 
positive and negative delays can be considered. So strictly 
speaking the cancellers are not truly real-time, the output is 
lagging in time by the length of the inserted delay.) So W is 
a vector with a complex element for each frequency channel. 
Figure [T| shows schematically how such a canceller can be 
implemented, hereinafter referred to as a mark one (MKla) 
canceller (the "a" is added since we will be modifying the 
filter later). 

Throughout this paper it is assumed that the filter weights 
are varying slowly enough that they are approximately con- 
stant on the time scales used to calculate them; typically less 
than a second. It is also assumed that the various signals are 
independent (i.e., Sm, Im, and 7V,„ are uncorrelated, and the 
receiver noise is independent for the different antennae, i.e., 
Nm and Nr are uncorrelated. The noise terms, N, will not be 
independent for receivers used in low frequency instruments 
such as LOFAR, where they will be dominated by partially 
correlated sky noise. Here we will only consider frequencies 
for which is dominated by uncorrelated noise, internal to 
the receivers.) When these assumptions hold, the power in a 
single frequency channel at the output of the canceller is 

P = {{Vm~WVr){Vrr,-WVr)*) 

= {\Vrr.?)-{v^w*v;)- ^3) 
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where * denotes a complex conjugate and (...) the expectation 
value. The explicit frequency dependency of the terms has 
been removed (it is assumed that the channels are independent 
so that Q can be applied to each frequency channel sepa- 
rately). If we seta? = a| = {\S\^), = (|iV,|2), 

and assume that the complex gains, delays, and weights are 
constant over the time average, so they can be taken outside 
the average, then we can write 



P 



w*g^g: 



|G„ 



2^2 



-j<l>n 



(4) 



where (j)rar represents the phase difference, 0„j — 0^- 

Equation (|4|i highlights a critical point. Minimising the out- 
put power reduces the combined power of the residual RFI 
and the inserted reference receiver noise. If there is refer- 
ence receive r noise the RFI will never be completely cancelled 
dWidrow & Stearns 1985; Barnbaum & Bradley 1998). This 
is because both the interference and reference receiver noise 
are being weighted. The receiver noise of the main antennae 
and the astronomy signals, however, are not affected by the 
choice of weights and will pass through the canceller freely 
(under assumptions of signal independence and zero reference 
gain towards astronomy sources). 

To perform the minimisation of P, one can differentiate it 
with respect to W and find the weights that set the deriva- 
tive equal to zero. The surface of P is a multidimensional 
(positive) quadratic surface that has a single minimum, so the 
weights that se t dP/dW = must give the unique global 
minimum ( Wid row & Stearns 1985; Barnbaum & Bradle"^ 
[1998). Alternatively, Wiener theory tells us that the optimal 
weights (known as the Wiener-Hopf solution), Wmkia, which 
minimise the output power are also the weights that set the 
cross-correlation between the canceller output and the refer- 
ence signal to zero C Widrow & Stearns. 1985. as indicated in 
figure [TJ: 
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So the weighting process takes the cross-correlation of the 
reference and main signals, which determines the correlated 
power and relative delay of the interference, and scales that 
by the auto-correlation of the reference signal. One can cal- 
culate these weights directly by calculating the correlations 
from short integrations, or they can be found adaptively by 
iteratively seeking and then tracking the weights that satisfy 
©. 

Ideally the output of the canceller would consist of the as- 
tronomy signal and receiver noise. However, as mentioned 
above, since there is always some receiver noise in the ref- 
erence signal, setting the correlation in (|5]l to zero can never 
remove all of the RFI. The RFI is played off against the ref- 
erence receiver noise. If INR^ is the interference-to-noise 
power ratio of the reference signal, G^aj/a"^ , the mean 
amount of residual output power (power in addition to the 
main receiver noise and the astronomy, Rmkia — Pmkia ~ 



Rmkl a 



(^2 n-2 

1 + INRr ' 



(7) 



It is clear from (|7| that as INRr approaches infinity (no 
reference receiver noise) the residual power goes to zero. If 
there is no RFI, will be zero and there will also be zero 
residual power (the filter turns off). When INRr is finite and 
non-zero the reference receiver noise term in the denomina- 
tor of ^ biases the weights and some residual power will 
remain. What might not be so clear from an inspection of 
(|7|l is the statement made earlier that this residual power is 
a combination of reference receiver thermal noise added dur- 
ing cancelling and residual RFI that was not excised. Another 
way to see this bias is to consider figure [T] Thermal noise 
from the reference receiver is present in both inputs to the 
weight generation process (in the Vr term from the reference 
antenna and the —WmkiaVr term from the filter output, c.f. 
equation|5]l. This will lead to a second non-zero-mean corre- 
lation product (the first being the RFI). Minimising the output 
power must be a trade-off between minimising the contributed 
reference receiver noise and the RFI, and as a result there will 
always be residual RFI. The total residual power given in dTJl 
can be divided into the inserted reference receiver noise resid- 

such that 



ual, R^rZL and the RFI residual, R^f,] 
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which are shown by lMitcheiill2004l to be 
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Thus R 



(rfi) 
mkl a 



R^Il\., I INRr. To shed some light on the 
meaning of the relations in (|9]) and dlOl l, we again interpret 
the process of output power minimisation as determining the 
weights that set the cross-correlation of the canceller output 
and the reference signal to zero. As INRr increases, the RFI 
becomes the dominant signal in the cross-correlation, and the 
filter must concentrate on reducing the RFI power As a result 
the proportion of the RFI that remains after cancelling de- 
creases faster than the injected noise power When INRr < 1 
noise starts to dominate and the filter will concentrate on re- 
ducing the contributed noise power When INRr goes to zero, 
the correlation is completely reference receiver noise, and the 
canceller turns itself off. When RFI dominates, most of the 
residual power is reference receiver noise, and when reference 
receiver noise dominates, most of the residual power is RFI. 
Also, since the reference signal is being scaled in an attempt 
to match its own RFI to the RFI in the main signal, the larger 
Gr is relative to Gm (for example using a reference antenna 
that is pointing directly at the interfering source), the smaller 
the scaling factor (weighting amplitude) and thus the amount 
of injected receiver noise. When INRr ^ 1 and Gr ^ Gm 
both of the residual terms drop off and extremely good results 
are achieved (Barnbaum & Bradley 1998; Bower 200 1|). 

2.3. MK2: Dual Reference Adaptive Cancellers 

To remove the biasing effect caused by the reference re- 
ceiver noise, a s econd reference re ceiver can be used, as 
was suggested in iBriggs et alj ( 120001) . The RFI in the main 
spectrum is still estimated using a weighted version of the 
spectral channels from the first reference, but now the cross- 
correlation of the second reference signal with the canceller 
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output is set to zero in order to find the weights. This MK2a 
canceller is shown in figure |2] 
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Fig. 2. — Schematic of a MK2a adaptive canceller, which uses two 
reference signals. Dashed lines indicate signal paths used to calcu- 
late the weights. 

It was mentioned in the previous section that the reference 
receiver noise is a component of both signals in Q. Zeroing 
this correlation results in minimum output power, but a non- 
zero RFI residual. Correlating the canceller output against a 
second reference with uncorrected receiver noise removes the 
bias. Since the RFI is the only correlated signal, zeroing the 
cross-correlation results in zeroing the RFI. However, there 
is no built-in mechanism to guard against the amount of ref- 
erence receiver noise contributed to the output. The amount 
added during cancelling must always be greater than the noise 
added by the MKla canceller, since power is no longer min- 
imised. 

The two references will be denoted ri and r2, so that 



Vr, = Nr, +G^,/eJ'^'-i,and 



(11) 



and we can again determine the weights that set the output- 
reference cross-correlation to zero: 



W„ik2a 



(V V*) 



(12) 



As with the MKla canceller the weights need to scale and 
phase shift reference signal Vr^ so that its RFI component 
matches that of signal Vm- Here, however, the independent 
reference signal r2 is used to give a true view of the scaling 
factor needed to match RFI levels. WmkZa is only constrained 
by the RFI, since the receiver noise terms in 
uncoiTelated, and the RFI signal is entirely replaced with a 
weighted version of thermal noise from reference signal Vr^ 
(assuming that the various signals are uncorrected and that 
the receivers are ideal and remain linear). While the noise in 
the weights will also increase the output power of the can- 
celler, it is much weaker noise since the weights are averaged 
over many samples, and is not considered further here. When 
Wmk2a from (fT2l l is substituted into the equation for output 
power, Pmk2a = {\Vm- Wmk2aVr^\'^) , the mean residual 
power is 
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as shown in lMitchell ( 120041) . Infinite attenuation of the RFI 
component has been achieved, but potentially a significant 
amount of system noise has been added during cancelling. 
Compare ^ and ([13]): As in (0, when INRr, > the at- 
tenuation of the RFI signal in Vm is very large. In this case 
however, any frequency channels with INRr, < 1 will end 
up with more unwanted power than they started with. The 
MK2a canceller does not turn itself off. The reference signal 
is boosted until its RFI matches the RFI in the main signal, 
regardless of the amount of receiver noise being added. Note 
that INRrg does not affect the output power, provided there 
is enough RFI power in Vr^ to keep the weights stable (see 
section[6]for a discussion of weight stability). 

If one is interested in the (interference-free) power spec- 
trum of Vm, the single reference antenna MKla canceller out- 
put will contain less residual power There are, however, ad- 
vantages to the MK2a canceller. If one is concerned with 
retrieving a structured signal from a voltage series, random 
noise in the signal may not pose too much of a problem, 
but a structured RFI residual may detract from signal recov- 
ery. More relevant in radio astronomy is the case where one 
is looking for a structure in the power spectrum. The RFI 
remaining after MKla cancelling will have features in the 
power spectrum, but if the (amplified and filtered) reference 
RFI, Gr-,(jj, is proportional to Gm<^j, equation ( fTSl ) says that 
the noise contributed by MK2a cancelling will have the same 
spectrum as the input reference receiver noise (apart from a 
constant scaling factor). It is also conceivable to remove the 
unwanted reference receiver noise from the auto-correlation 
of either canceller's output. This is discussed in the next sec- 
tion and can result in a MK2 canceller superior to the MKl 
canceller for some applications. 

2.4. Suppressing Added Reference Noise 

We now describe a double canceller setup that can be used 
to suppress the added reference receiver noise in the output 
astronomy power spectrum. If the main signal is duplicated 
before cancelling so there is an identical copy, the copy can 
be passed through a second canceller that uses different refer- 
ence signals so that the noise added will be uncorrected with 
the noise added to the first. When the two filtered copies are 
cross-coiTelated, the main antenna receiver noise and astron- 
omy will correlate as if the original signal has simply been 
auto-correlated, while the added reference noise power will 
average away with the radiometric factor {\/ Ai/Tint)- 

An important point to note about the MK2a canceller is that 
all of the noise added during cancelling is from the first refer- 
ence receiver (see equationfl3]). The second reference is only 
used in setting the complex weights. So the second canceller 
can be made by interchanging the references. This setup will 
be called the MK2b cance ller and is shown in figure[3] As dis- 
cussed in lMitchell (120041) . if the INR of the references are the 
same, the mean residual output power from the MK2b can- 
celler is 



R 



R 



rnk2a 



mk2b — 



(14) 



which averages towards zero as the integration length is in- 
creased. 

Similarly, if two references are used to create two indepen- 
dent MKla cancellers for the two main signal copies, this 
gives the MKlb canceller shown in figure U and the mean 
residual output power becomes 
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Fig. 3. — MK2b canceller system where the auto-correlation is 
replaced by the cross-correlation of two independently processed 
copies of the main signal. 
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Fig. 4. — MKlb canceller system where the auto-correlation is 
replaced by the cross-correlation of two independently processed 
copies of the main signal. 

Equations ST4i and ( fTSl ) show that while RmkBb should keep 
integrating towards zero, Rmkib has a definite limit due to the 
RFI signal that remains after cancelling. However, one must 
be aware that in situations where the reference INR is very 
small the MK2 canceller does not turn itself off, and there 
are practical implementation issues that need to be addressed 
(essentially, one may need to force the canceller to turn off). 
This is discussed in section |6] 

3. THEORETICAL RESULTS 

In this section the residual power equations given through- 
out section |2] are demonstrated. Note that while the plots for 
the dual-reference cancellers show the residual power as the 
reference interference-to-noise power ratio approaches arbi- 
trarily close to zero, the algorithms in practice become unsta- 
ble and need to be turned off. This point will be reiterated 
where appropriate in the discussion below. For the theory we 
have set (7| — aj^ = 0, so that all of the output power dis- 
played in this section is a combination of residual RFI and any 
reference receiver noise added during cancelling. 

Figure |5] displays the proportion of residual power in the 
output signal after adaptive cancelUng with MKlb and MK2b 



cancellers. The plot shows that the added reference receiver 
noise averages away as the integration length is increased. If 
single canceller systems were being considered, then since 
any noise added during cancelling is sent to an auto-correlator, 
the output power would remain constant (i.e., remain at the 
levels shown on the left hand side of figure|5]). It is clear that 
the MKlb canceller hits a limit when it reaches the residual 
RFI, but that the MK2b does not. 
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Fig. 5. — Theoretical residual power (normalised by the main sig- 
nal's input RFI power) for various INRr values as a function of the 
number of samples used to average the noise down in the canceller 
output (AuTint). The three sets of solid and dashed lines are for 
INRr values of 0.1 (top), 10 (middle), and 1000 (bottom). Solid 
lines indicate MKlb cancelling, dashed lines MK2b cancelling. The 
lines flatten when the residual RFI power level is reached. 

As INRr decreases the normalised output power of the 
MKlb canceller levels off at 1, so there is no cancelling tak- 
ing place. On the other hand the MK2 canceller continues to 
insert more and more reference receiver noise in an attempt to 
match the reference RFI to the main signal RFI. Even though 
the MK2b canceller always has the larger total residual power, 
it is entirely zero-mean noise and averages out with the radio- 
metric factor. Again the reader should note that for low INRr 
values the MK2 canceller can become unstable and requires 
an additional mechanism to turn off. 

Figure |6] shows contours of constant (normalised) output 
power as a function of INRr and the number of samples, 
AvTint- Figures [6^ andjSJ) represent MKla and MK2a can- 
celling respectively, and figures |6}; and|6}i represent MKlb 
and MK2b cancelling respectively. The amount of residual 
power in dB is indicated by the grey scale and runs from -80 
to 20 dB. It is clear from figure |6}; that a constant RFI resid- 
ual remains for all INRr values after MKlb cancelling. The 
dashed line indicates the approximate line where the added 
reference receiver noise power has averaged down to expose 
the non-zero residual RFI power level. 

It is clear from figures [6^ and |6}; that an output power 
plateau is reached as INRr becomes small for the MKl can- 
cellers. The residual power of the plateau is dB, and indi- 
cates that the canceller has turned off. In contrast, the MK2a 
canceller ^jp) does not turn off and results in more output 
power than input RFI power for low INRr levels. However, 
since the residual is entirely noise it averages down in a MK2b 
canceller, where two independent filters are used (|6}l). This is 
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Fig. 6. — Residual power contours as a function of INR,. and lS.VTint. Greyscale indicates the normalised output power, Rmm/{\Gm\'^cr'j), 
expressed in dB. Figures a and b show MKla and MK2a cancelling respectively. Figures c and d show MKlb and MK2b canceller output 
respectively. 



highlighted further in the next section. 

4. MICROWAVE RFI AT THE ATCA 

We now demonstrate adaptive cancellation of real RFI im- 
pinging on the Australia Telescope Compact Array (ATCA). 
The RFI is a point-to-point microwave (MW) television link 
transmitted from a TV tower on a nearby mountain at 1503 
MHz. The reference antennae were two orthogonal linearly 
polarised receivers on a small reference horn poi nted in the 
direction of the MW transmitter, as described in iBell et afl 
(|2001)Q A linearly polarised receiver on a regular ATCA an- 
tenna, pointing at the sky and receiving the microwave link 
interference through the antenna side-lobes, was used to col- 
lect the main signal. The RFI is polarised, and as long as all 
three receivers are at least partially polarised in the same sense 
as the RFI the cancellers will work correctly (assuming that 
the polarisation cross-talk between the reference receivers is 
negligible so the receiver noise is independent). The received 
voltages were filtered in a 4MHz band centred at 1503 MHz, 
downconverted, and sampled with 4-bit precision. 

Each of the cancellation techniques discussed has been 
applied to the M W data using MATLAB (see the 
ImATLAB User's Gu ide 1998). All of the spectra shown were 
generated using 1024-point FFTs. Figure [T] shows the unfil- 
tered and MKla filtered power spectra of the main astronomy 
voltage series when large reference INRs were available. Two 

Dataset srtca02. 



synthetic astronomy signals have been added to the astron- 
omy voltage series at 1503.0 MHz and 1503.2 MHz. What 
we want is to remove the RFI peak (the 1-2 MHz wide peak 
centred at 1503 MHz) from the main signal spectrum, leaving 
the broad, main antenna bandpass, and the astronomy. This is 
indeed what is seen, and similar results are achieved for all of 
the techniques discussed in this paper. Since all of the tech- 
niques behave excellently when the reference INR is large, 
and any contributed reference receiver noise or residual RFI 
is much smaller than the main signals receiver noise level, it 
is difficult to compare them. 

Figure|8]displays the power spectra of two reference signals 
that have had Gaussian random noise injected into their sam- 
pled voltage series, to set INRr « J at the centre of the MW 
band and zero at the edges. The reason for adding the fake 
receiver noise was to lower the INR and accentuate any resid- 
ual RFI, as detailed in Q, dTOl l. and ( flJl l. Plots of the power 
spectra of the main signal (the one containing the astronomy), 
before and after the different cancellation techniques, are dis- 
played in figures l9lthrough[T2l 

Figure |9] shows unfiltered and MKl filtered power spectra 
of the main astronomy voltage series. The two synthetic as- 
tronomy peaks have not been affected to any measurable level. 
The amount of residual power after MKla cancelling is as ex- 
pected from (|7]l. About twice as much power has been can- 
celled at the centre of the RFI peak by the MKlb canceller 
(the added receiver noise in this case has a zero-mean cross- 
correlation), which is expected since INRr ~ 1 ■ Away from 
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Fig. 7. — Power spectra of main signal Vm before and after MKl cancelling 
with INRr ^ 1. Two simulated cosmic signals were added at 1503.0 and 
1503.2 MHz. 



the centre of the peak, INR,. becomes smaller and the propor- 
tion of residual RFI increases (which is why the residual RFI 
peak is flatter than the initial peak). No matter how long the 
integration is run this residual power will not decrease any 
further. That is, the canceller is operating in the horizontal 
region shown in figure |5] 

Before and after spectra for the MK2 cancellers are shown 
in figure[TO] As for the MKl cancellers, the two synthetic as- 
tronomy signals have been left unaffected. It is clear from 
figure [foji that when INRr < 1 the MK2a output resid- 
ual power is greater than the input RFI power (except when 
INRr goes to zero, where the canceller was turned off, see 
section |6]l. This is clearly unsuitable and makes the situation 
worse. In contrast, the MK2b canceller (figure [TOb) has re- 
moved the RFI peak down to the primary antenna's receiver 
noise level, indicating that the MK2a residual power was com- 
pletely added noise and not residual RFI. The removal of the 
RFI was extremely successful, with only a slight increase in 
zero-mean noise at the primary receiver noise floor, as seen in 

figure nn 

Each spectral channel in the examples in this section 
was averaged for 15872 samples. This means that for the 
MKlb and MK2b cancellers the contributed reference re- 
ceiver noise power will have averaged down by about a factor 
of 715872 w 126, or about 20 dB. Since the MK2b can- 
celler output has essentially the same noise spectrum as the 
MK2a canceller output, but averaged down, there has been 
about a 20 dB reduction in unwanted power after MK2b can- 
celling at the centre of the RFI peak (see figure [TOl l. This 
zero-mean noise would have continued to average down had 
we continued to integrate, and more attenuation would have 
been achieved. The MKlb canceller has only achieved around 
6 dB of attenuation, and this will not increase with a longer 
integration. One should keep in mind, however, that we have 
inserted a substantial amount of noise into the reference volt- 
ages, which has also greatly lowered the achieved attenuation. 

An important final point to note is that, with the exception 



of the edges of the RFI peak (where there are features associ- 
ated with turning the filter off), figure [TOh shows that the RFI 
is spread with essentially constant power over the RFI con- 
taminated channels. This is because INRr, which we recall is 
given by G^ctj / , is approximately proportional to G^^^aj, 
so that the spectral shape of the RFI cancels out of the residual 
power given in ( fTsl i. 

5. POST-CORRELATION CANCELLERS 

From an implementation point of view it is important to 
realise that in astronomy we are not usually seeking fast time- 
variable information such as modulation. In most astronom- 
ical applications the aim is to measure signal statistics since 
they are related to quantities such as cosmic flux density and 
visibilities as measured by arrays. These applications are gen- 
erally either finding the auto-correlation of signals from a sin- 
gle antenna (to measure the power spectrum of the astronomy 
signal), or the cross-correlation of signals from more than one 
antennae (to measure the spatial coherence of the astronomy 
signal). Se e Thoms on et al.l (11986) and various chapters of 
iTaylor et al.i (.1999) for an overview. In this paper we have 
concentrated on the power spectra given by auto-correlations, 
however all of the techniques discussed can be generalised to 
work on cross-correlations. 

The upshot of only requiring signal statistics is that if the 
canceller weights are not changing appreciably over the 100 
millisecond or so time interval that the statistics are measured 
over, then the algorithms can be applied to the statistics rather 
than each voltage series. This means that they are applied at a 
rate of Hz to kHz rather that MHz to GHz. It also means that 
if the new astronomical cross-correlators that are coming on- 
line for new and existing facilities have a few extra inputs for 
reference antennae, then no new filters need to be added to the 
signal paths, and the cancellation can be performed after the 
observation as part of post-processing. For a comparison of 
voltage cancellation and statistics cancellation for RFI signals 
that require filters with weights that are changing appreciably 
during an integration, see Mitchell & Robertson (2005). 

The standard statistics canceller used in radio astronomy is 
a post-correlation version of the MK2b canceller discussed in 
section 1273] If we use Pjk to denote the correlation between 
signals from antennae j and k, the quantity that is subtracted 
from the auto-correlation of antenna m is determined from 
amplitude and phase closure relations (fBriggs et alJl2000h : 



pc2 



P P* 



\Gr. 



|2^2 

CTr. 



(16) 



Since the reference signal does not contain any information 
about the astronomy signal, it will not be present in any of 
the quantities in ( fT6b . The denominator removes the refer- 
ence signal RFI's phase and amplitude information from the 
numerator, leaving information about the RFI in the signal 
from antenna m, and zero-mean noise, since the expectation 
operators are not infinite in extent. When Cpc2 is subtracted 



from main si gnal power P„ 
lMitchellll2004i) 



Rpc2 — 



1, the mean residual power is (see 



^INRrJNRr.^/Aunnt' 



(17) 



as in (fl4] i. Rpc2 is shown in figure [12] As long as the gain 
and geometric delay stay essentially constant over the time 
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Fig. 8. — Power spectra and estimated receiver noise of the reference signals. Most of the bandpass noise was artificially added to lower 
the INR. The peak at 1502.7 MHz is an artifact of the receiving system that leaked into the reference voltages, (a) Reference signal K-i- (b) 
Reference signal Vrj • 
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Fig. 9. — Power spectra of main signal Vm before and after MKl cancelling. Two simulated cosmic signals were added at 1503.0 and 1503.2 
MHz. (a) MKla canceller, (b) MKlb canceller. 



average, the "pre" and "post" correlation techniques are very 
similar; Rpc2 ~ Rmk2b- 

While in the last few sections an extreme case has be 
demonstrated (that of a low reference INR), it does highlight 
major differences between the different cancellers. The main 
differences being that infinite RFI removal is theoretically 
possible for cancellers that use two reference signals, but not 
for single reference cancellers. It is also clear that when the 
weights are slowly varying, the post-correlation canceller is 
equivalent to the dual-reference MK2b canceller. 

Another way in which the MK2b and post-correlation can- 



celler techniques can differ is in their instabilities at low INR 
levels. Equations ( fT2] l and ( fT6b have denominators that be- 
come zero-mean noise when there is no correlated RFI signal, 
which can lead to numerical errors. The coefficients of the 
single-reference techniques go to zero when the RFI becomes 
weak and the cancellers automatically turn off. These issues 
are discussed next. 

6. INSTABILITIES IN THE DUAL REFERENCE ALGORITHMS 

Equations ( fT2] i and ( fTSI l show that the dual-reference can- 
cellers can have stability problems when INRr <C J . In fre- 
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Fig. 10. — Power spectra of main signal Vm before and after MK2 cancelling. Two simulated cosmic signals were added at 1503.0 and 1503.2 
MHz. (a) MK2a canceller, (b) MK2b canceller. 
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Fig. 11. — MK2b canceller spectrum from figure [ToJ), with a INB^r > 
1 MK2b spectrum overlaid. Very little noise is added when INRr ^ 1, 
demonstrating the increase for the INRr < 1 case. 

quency channels where the correlated interference is zero or 
very small, (fT2] i and ( fTSI l are noise dominated and can result 
in a division by zero (or very close to zero). This cannot oc- 
cur in the MKl cancellers since (|6]l goes to zero as / — > 
and they turn off. A modified post-correlation canceller has 
been suggested by Briggs et al. (2000) in which an extra term 
is added to the denominator of (fTSI l: 




r" 



P P* P 

rnr\ rnr2 r\r2 



*2 



P* P 



(18) 
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Fig. 12. — Power spectra of main signal Vm before and after post- 
correlation cancelling. Two simulated cosmic signals were added at 1503.0 
and 1503.2 MHz. There are slight artifacts at the edge of the peak where 
INRr — ^ and the canceller is turning off, as described in section|6] 

reduces to (fT6] l. 

Since introduces a small bias in a similar way to the 
MKl cancellers, a relation equivalent to dTji, but with a much 
smaller bias, can be derived: 



where * = {(K^K^) (K^NQ /Ai^tm)'^\ which is an 
estimate of the noise power in the reference cross-correlation, 
and the prime indicates that Cpc2 has been approximated. The 
extra term ^I*^ stops the zero-mean fluctuations in -P^ir, from 
going too close to zero, while for large INRr equation ( fTSl l 



'pc2 



l + INRiINRsAvnnt 



(19) 



This will have both a noise and a RFI component, as in 
(|7l), however, since the canceller is now biased like a MKl a 
canceller the added receiver noise term will not average away. 
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R'pc2 reduces as INR^ multiplied by the number of samples 
in the time average. 

There is a similar problem for the pre-correlation MK2 
adaptive cancellers. In the lag domain, the division in 
([12] ) becomes a multiplication by the inverse of a matrix 
with columns containing offset copies of th e ri-r2 cross- 
correlation function (Wid row & Stearnslll985l) . Divisions by 
zero in the frequency domain due to interference-free fre- 
quency channels in P^^r-, ^re manifest in the reference lag 
matrix as singular values. One method of dealing with this 
is to use singular value decomposition to decompose the ma- 
trix into two orthonormal tria ngular mat rices and one diago- 
nal matrix (see § 2.9 of Press et al.lll986 ). Singular (or near- 
singular) values can be selected when the relevant diagonal 
matrix elements are less than a chosen threshold, such as ^ . 
The singular parts of the matrix contain no information about 
the correlated signal and are removed from the decomposition 
matrices. The inverse matrix can then be reconstructed from 
the remaining parts of the three decomposition matrices, and 
it will not function in the RFl-free parts of the spectrum. 

7. SUMMARY 

Interference cancelling using a single reference antenna can 
give excellent results when the reference signal interference- 
to-noise ratio is large, and there is more gain towards the in- 
terfering signal for the reference antenna than for the astron- 
omy antennae. However, receiver noise in the reference signal 
means that a fraction of the interference will always remain 
after cancelling. A second reference signal can be used to re- 
move the noise bias and give infinite interference attenuation, 
but a larger amount of reference receiver noise is added during 



cancelling. For pre-correlation systems, a dual canceller setup 
can be used to average the (zero-mean) receiver noise away, 
a process that comes automatically with post-correlation can- 
cellers. A breakdown of the main properties for the different 
mitigation techniques is given in table m 

It is important to note that even though the single-reference 
cancellers leave residual RFI, the residual may be extremely 
small and well below the RMS noise. This occurs when the 
interference-to-noise ratios of the reference signals are very 
large, and the use of two references (in pre-correlation sys- 
tems) might just add complexity to the system with little or 
no benefit. However, if maximum sensitivity is required, one 
should be aware that they will eventually reach a non-zero 
residual signal. 

Using two reference signals to remove the reference re- 
ceiver noise bias removes the inherent stability of the algo- 
rithms in situations where some or all of the frequency chan- 
nels are interference-free. Although there are applications in 
which the passband will always be entirely filled with RFI 
(such as observations in the GPS LI and L2 bands), many in- 
terfering signals will only take up a part of the band. In these 
cases the algorithms need a mechanism to turn themselves off 
in the vacant frequency channels. 

We are grateful to Dr Mike J. Kesteven and Professor 
Lawrence E. Cram for discussions and comments on this pa- 
per. The Australia Telescope Compact Array is part of the 
Australia Telescope which is funded by the Commonwealth 
of Australia for operation as a National Facility managed by 
CSIRO. 
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Fig. no. 


Ref. 
Ants. 


Filters 


Configuration 


Behaviour 


Eqn. no. 


MKla 


m 

LU 


1 


1 


One reference antenna feeds both fil- 
ter and weight generator. 


Minimises residual power, but in- 
cludes some RFI. Turns off gracefully 
as INRr 0. 


I7II9II10I 


MK2a 


|2| 


2 


1 


Independent reference antennae for 
filter and weight generator. 


Totally cancels RFI, but adds ex- 
tra noise. Must be controlled as 
INR, ¥ 0. 




MKlb 


|4| 


2 


2 


Two MKla cancellers, operating on 
copies of main signal. 


Sensitivity decrease, and non-zero 
RFI floor. 


li5j 


MK2b 


13J 


2 


2 


Two MK2a cancellers, operating on 
copies of main signal. 


Sensitivity decrease, but no non-zero 
RFI floor. 


iMl 


PC2 


not 
shown 


2 


1 


Two reference coirelations with each 
astronomy signal. 


Sensitivity decrease, but no non-zero 
RFI floor. Must be controlled as 
INRr 0. 
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TABLE 1 

Summary of the main properties for various interference mitigation techniques 



